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TITLE: 

TYPE: 

NUMBER OF WORDS: 

TEMPORARY STORAGE: 

ACCURACY: 

DURATION: 

DESCRIPTION: 



NOTE 1: 
N0TE:2: 



Error Function (DOI or SADOl) 

Closed, standard entry 

22 

0, 1, 2 

-39 / \ 

Maximum error + 10 x 2 (see below; 

10 milliseconds 

This routine computes the function 



erf (x) = 



2 f* 



e dt 



for |x] < 1 by use of a polynomial expansion. The 

routine is to be entered with the argument in the 

accumulator and the link address in the right address 

position in Q. When the link is obeyed, the computed 

function will be in the accumulator and quotient 

registers <> 

An incorrect result is given when x = -1* 

The rounding of the coefficients in the power series 

is done in sueh a way that the error in the computed 

values will be found to be considerably smaller than 

the maximum given above. The errors for a large range 

of experimental values were found to be less than 

-12 

5 x 10 ^. 



nj 



DAT E December 12, i960 
PROGRAMMED BY- 
APPROVED BY 




LOCATION 


OBDER 


NOTES PAGE 1 


v 13 




00K(V13) 









4o F 
K5 13L 






1 


42 10L 








90:. -f 


Plant link 




2 


7J^F 








40 IF 


x at 1 




3 


L5 L 








42 5L 


Plant coefficient address 




4 


23 5L 








50 2F 


Enter loop 




5 
6 


7J 1? 
L4 F 
40 2F 
F5 5L 


8 

Form Z ^ C x 
n=0 n 




7 


42 5L 








LO 11L 


End test 




8 


32 4L 
50 2F 






9 


75 F 








66 12L 


Form x(E C n x 2n ) -7 ^2 




10 


S5 F 








22 F 


Exit via link 




ll 


U IF 








L4 22L 


End constant 




12 


40 F 








00 38 6226 925 453J 


V^/2 




13 


00 F 








00 93 8982 J 


Cg = .000000938982 




111 


NO F 








00 4999 8767 1076J 


c = -.000012328924 




15 


00 F 








00 1059 6568 2J 


C£ = .000105965682 




16 


NO F - 








00 4992 4292 5533J 


c = -,000757074467 





LOCATION 



ORDER 



NOTES 



PAGE 2 



V 13 



17 
18 

19 
20 

21 



00 F 

oo 4629 ^589 26J 

NO F 

00 V761 9050 8999J 

00 F 

00 9999 9996 832J-; 

NO F 

00 1666 6666 6785J, 

80 F 

00 F 



C = o 004629458926 



c = -.023809491002 



C 2 = .0999999968.32 



G 1 = -o333333333215 



C =1 






